Kondo-lattice model: Application to the temperature-dependent electronic structure of EuO(lOO) 

films 
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We present calculations for the temperature-dependent electronic structure and magnetic properties of thin 
ferromagnetic EuO films. The treatment is based on a combination of a multiband-Kondo lattice model with 
first-principles TB-LMTO band structure calculations. The method avoids the problem of double-counting 
of relevant interactions and takes into account the correct symmetry of the atomic orbitals. We discuss the 
temperature-dependent electronic structures of EuO(lOO) films in terms of quasiparticle densities of states and 
quasiparticle band structures. The Curie temperature Tc of the EuO films turns out to be strongly thickness- 
dependent, starting from a very low value (~ 15K) for the monolayer and reaching the bulk value at about 25 
layers. 



i 

C 

o 
o 



(N 
> 

in 

(N 
O 



c3 



O 

o 



X 



I 



PACS numbers: 75.50.Pp,73.20.At,75.70.Ak,71.15.-m 



I. INTRODUCTION 

A highly vivid field of research in current solid state 
physics aims at the intercorrelation between electronic and 
magnetic properties of materials. Special highlights in this 
respect include phenomena like the colossal magnetoresis- 
tance (CMR) effect observed in the manganites perowskites 
Ti- x D x Mn03 where T is a trivalent lanthanide ion such 
La and D is a divalent alkaline-earth (e.g. Ca, Ba, Sr) ionl 
The interest is fueled by the prospect for applications in the 
fields of storage media and sensor technology. Interesting 
from an even more fundamental point of view is the intercor- 
relation of electronic and magnetic properties in the so-called 
iocal-moment' magnets, protolypically realized by the metal- 
lic rare earths Gd, Tb ancLDyB and the insulating rare-earth 
compounds EuO and EuStl. Sometimes these materials are 
also denoted as '4f systems' because a great part of their prop- 
erties is due to the existence of a partially filled and highly lo- 
calized 4f-shell belonging to the rare-earth ion. Strong screen- 
ing by other closed electron shells keeps the Hund's rules 
valid. According to these the 4f electrons couple to a finite 
moment which is strictly localized at the sites of the rare-earth 
ion (Eu 2+ , Gd 3+ , ...). While the purely magnetic properties 
of these materials are due to their localized 'magnetic' 4f elec- 
trons the conductivity properties are determined by the quasi- 
free electrons in rather broad (5d, 6s) conduction bands. 

Many characteristic features of 4f systems follow from an 
intimate correlation between localized and itinerant electrons. 
This correlation leads, e.g., to a striking temperature depen- 
dence of the conduction band states induced by the magneti- 
zation state of the local moment system. First evidence for 
this has been found in 1964 by the discovery of the 'red shift 
effect' in EuO by Busch and WachteiO. The effect refers to 
a strong shift of the lower band edge to lower energies upon 
temperature cooling below To A further striking effect be- 
ing due to the mentioned induced temperature dependence is 
a metal-insulator transition observed in Eu-rich EuO-J. The 
Eu-richness manifests itself in (2+)-charged oxygen vacan- 
cies. One of the two excesses electrons is thought to be tightly 



bound to the vacancy while, because of the Coulomb repul- 
sion, the other occupies an impurity level fairly close to the 
lower band edge. When cooling below Tc the band edge 
crosses the impurity level ('red shift') freeing therewith the 
excess electrons. A conductivity jump of up to 14 orders of 
magnitude has been observed. 

Further effects result from the interaction of the band elec- 
tron with collective excitations of the moment system. One 
of these is the creation of a new quasiparticle which is called 
'magnetic polaron'. It can be classified as a propagating elec- 
tron dressed by a virtual cloud of excited magnons. The for- 
mation of magnetic polarons strongly influences the electronic 
structure of respective materialsEDQ. 

The influence between localized and itinerant electrons is 
of course of mutual character. That means that a finite band 
occupation has observable consequences for the moment sys- 
tem, too. The ferromagnetic order in Gd, e.g., can be ex- 
plained only by indirect exchange coupling (RKKY) of the 
localized moments (spins) mediated by a spin polarization of 
the a priori non-magnetic conduction electrons. This RKKY- 
mechanism can excellently be monitored by inspecting the al- 
loy Eu\- x Gd x S . Replacing the Eu 2+ -\on in the ferromag- 
netic insulator EuS by a Gd 3+ leads in a definite manner to a 
population of the conduction band without diluting the local 
moment system. The latter holds because Eu 2+ and Gd 3+ 
both have a half-filled 4f-shell S = |. That allows for a di- 
rect study of the carrier concentration dependence of the Curie 
temperature T c = T c (n)EB. 

It is commonly accepted that the so called s-f modeOLil 
provides a good starting point for an at least qualitative un- 
derstanding of the above properties of typical 'local moment' 
materials. Within the model they are traced back to an in- 
traatomic, interband exchange interaction between localized 
and extended states. The same model has been extensively 
applied in the recent past to the CMR-materials. In the man- 
ganites the three almost localized 3d t2g -electrons couple to an 
S = | -spin that is exchange-coupled to the somewhat more 
mobile 3d eg -electrons. In this context the s-f model is referred 
to as the ferromagnetic Kondo-lattice model (FKLM)I131jo. 
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The additive 'ferromagnetic' points to an exchange coupling 
which favors the parallel alignment of itinerant and localized 
spin. Without this additive the 'normal' Kondo-lattice model 
(KLM) is meant that differs from the FKLM only by a sign 
change of the coupling constant favoring therewith an antipar- 
allel alignment. The latter version is considered a candidate 
for the extraordinary heavy-fermion physicslid In this paper 
we exclusively refer to the FKLM, which covers the above de- 
scribed 4f-systems and the manganites. While the 4f-systems 
surely belong to the weak-coupling region (exchange coupling 
much smaller than the bandwidth) the manganites are to be as- 
cribed to the strong coupling limit. 

A central aspect of our study which we are going to present 
in this paper is that of reduced dimensionality and its interplay 
with typical correlation effects. Magnetic phase transitions 
and electronic structures at surfaces and in thin films have at- 
tracted a lot of recent research work. One of the most remark- 
able magnetic phenomena observed in 4f-systems is the ex- 
istence of a ferromagnetic surface at temperatures where the 
bulk material is already paramagnetic. This effect was, first 
documented for the Gd(0001) surface by Weller et al.O and 
since then has been measured by several groups with different 
experimental techniques. Depending on the used techniques 
and maybe on the cleanness of the applied surface Curie- 
temperature enhancements ATc — Tc (surface!-,— Tc(bulk) 
in between 17K and 60K have been foundHrl 1 '!. However, 
the possibility of an enhanced surface Curie-temperature re- 
mains a matter of controversy. Donath et alJij using spin- 
resolved photoemission did not find any indication for ATc > 
at Gd(0001) surfaces. By careful analysis applying spin- 
polarized secondary-electron emission spectroscopy to sur- 
face magnetization measurements and magneto-optic Kerr ef=. 
feet to bulk magnetization measurements Arnold and PappasEj 
come to the same conclusion. With respect to a possible ATc 
at Gd (also Tb) surfaces the existence of a Gd(0001) sur- 
face state seems to nlay a crucial role. Predicted by band 
structure calcuJaJianscJ and simultaneously measured with 
photoemissionEjoEfl in particular its temperature dependent 
induced exchange splitting has attracted attentions. Is it 
'Stoner-like' collapsing for T — > Tq or does the splitting 
persist in the paramagnetic phase performing the demagne- 
tization by_-a_redistribution of spectral weight ('spin-mixing 
behavior ')EiE3? It is an open question whether and how a 
surface state can cause an enhanced surface-Curie tempera- 
ture. To get insight a calculation of the temperature-dependent 
electronic structure of the real (!) Gd-film would be required. 
For providing the clues to such a treatment of local-moment 
metals we present in this paper the methodical approach for 
dealing with local-moment semiconductor films, which are in 
addition interesting in themselves. The semiconducting 4f- 
systems EuO, EuS may be considered as 'low-density limits' 
of the metallic Gd. They have the same half-filled 4f-shell as 
Gd producing localized S = |-spins. They differ by the car- 
rier concentration n (n = in EuO, EuS; n ^ in Gd). We 
present in this paper the description of a real EuO-film of vari- 
able thickness. We are interested in the layer-dependent mag- 
netic and electronic properties. For this purpose we combine 
a many-body treatment of the FKLM with a 'first-principle' 



electronic structure calculation. For the many-body part we 
use a 'moment-conserving decoupling procedure' (MCDA) 
the idea of which is developed in ref. 6 for the single-band 
bulk-FKLM. The reformulation for the-special situation of 
a film geometry has been performed inll We present here 
the generalization to a multi-band-FKLM film. As to the 
real EuO (100) film we combine the model calculation with 
a 'tight-binding linear muffin tin orbital' (TB-LMTO) band 
structure calculation to get temperature- and layer-dependent 
magnetizations, quasiparticle densities of states (Q-DOS) and 
quasiparticle band structures (Q-BS). We find for thicker films 
(> 20 monolayers) a surface state. Being of course unoc- 
cupied it nevertheless exhibits a temperature-dependent ex- 
change splitting induced by exchange coupling to the mag- 
netic 4f-moment system. We speculate that this surface state 
can be considered an analogous to the Gd-surface state men- 
tioned above. Furthermore, we find some evidence that the 
temperature behavior of the EuO-surface state together with 
the red shift of the lower conduction band edge may close the 
4f-5d gap in the j-spectrum giving rise to a surface-insulator- 
halfmetal-transition. 



H. MULTIBAND-KONDO LATTICE MODEL 
A. Model-Hamiltonian 

We investigate a film consisting of n equivalent layers par- 
allel to the surface of the film. Each lattice site of the film 
is indicated by a Greek letter a, (3, j, . . . denoting the layer 
index and a Latin letter i, j, k, . . . numbering the sites within 
a given layer. Each layer possesses two-dimensional trans- 
lational symmetry. Accordingly, the thermodynamic average 
of any site dependent operator Ai a depends only on the layer 
index a: 

{A ia ) = (A a ) . (1) 

The complete d-f model Hamiltonian for a real system with 
multiple conduction bands, 

H = TLd+'Hf +Hdf, (2) 

consists of three parts. The first 

7~L<1 = ^ija}} C tamcr C jPm'a (3) 

ijafltr mm' 

contains the conduction band structure of e.g. EuO. The in- 
dices 77i and m! denote the different conduction bands. cf ama 
and Ci ama are, respectively the creation and annihilation oper- 
ators of a spin-er electron from the m'th subband at the lattice 
site R Jct . T£-™g are the hopping integrals which later have to 
be obtained from an LDA calculation. 

Each lattice site Ri a is occupied by a localized magnetic 
moment, represented by a spin operator Sj a . In the case 
of Eu and Gd the 4f shell is exactly half-filled and, be- 
cause of Hund's rules, has its maximal magnetic moment of 
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7 fiB- These localized moments are described by an extended 
Heisenberg Hamiltonian 



fit 



ija/3 



(4) 



The first term is the original Heisenberg Hamiltonian with 
the exchange integrals Jff ■ The second term introduces a 
symmetry-breaking single-ion anisotropy, which is necessary 
for film geometri es to obtain a collective magnetic order at 
finite temperaturesc^EZI. Here, one has typically Dq <C J"/- 

The distinguishing feature of the single-band Kondo-lattice 
model (s-f or d-f model) is an intraatomic exchange between 
the conduction electrons and the localized f-spins. The form 
of the respective Hamiltonian in the case of multiple conduc- 
tion bands can be derived from the general form of the on-site 
Coulomb interaction between electrons of different subbands. 



L\. . .Zj4 <j(j' 



'^2u Ll ... L4 c+ ia c+ 2ljl c L3a ,c L4a (5) 



We drop for the moment the site index i because we want to 
keep the local character of the s-f interaction. L1...L4 de- 
note the different bands, and Ul 1 ...l 1 are the Coulomb matrix 
elements. Restricting the electron scattering processes caused 
by the Coulomb interaction to two involved subbands we get 
instead of (||): 

Hl = \ Yj Y^{ UL - L ' C tv C l'v' C L'°' C L° + 



L,L' 



Jl,L' c Lcr C ia' C L'cr' c L'a] 



(6) 



ant electrons: 



H c /- ^ 



/ , [Umf c mcr c fcri c fcr' c ma + 
7n,f,<r,<r' 

+ 



1 



2^mf C ma C ma' C f^' C f cr + 



2 J fm C fa C fa' C 



ma' wner J 



(8) 



The last two terms do not contribute in case of the 4f systems 
(EuO,Gd) since the (Eu 2+ , Gd 3+ )-4f shell has its maximum 
spin of S = |. All the seven 4f electrons have to occupy 
different subbands and none of the 7 sublevels will be doubly 
occupied. Introducing electron spin operators, 



a 

a' 



(9) 



one arrives at 



mf 

+ y^XUmf - ^Jmf) n m Uf (10) 
mf 



with n 



m(f) 



In the case of EuO, we 



are interested in, the last term does not contribute because 
of n m = 0. Even for Gd (n m ^ 0) it certainly can be ne- 
glected since the number of electrons is fixed, n f thus being 
only a c-number. By defining the spin operator S of the local 
f-moment, 



Thinking of 4f systems such as EuO and Gd the band indices 
L and L' can be attributed either to a flat 4f band (L — > /) or 
to a broad (5d,6s) conduction band (L — > to). In an obvious 
manner we can then split the Coulomb interaction into three 
different parts 



Hi = HT + H f / 



(7) 
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and by assuming the interband exchange J m f to be indepen- 
dent on the band indices m and f, 



J rnf = - ^SJ 



(12) 



depending on whether both interacting particles stem from 
a conduction band, Jif, or both from a flat band,7Yj , or 
one from a flat band the other from a conduction band,7i^. 
The first term refers to electron correlations in the conduc- 
tion bands, which by definition is neglected in the KLM. In 
the case of EuO this simplification becomes exact because the 
conduction bands are unoccupied (n = 0). On the other hand, 
for the CMR-materials the neglect of H.f appears question- 
able. Some authors therefore extend the KLM to the 'corue.- 
lated' KLM by regarding TUf as the Hubbard-interactionEl 
The second term,7i^, is already contained in the description 
of the localized moments via the Heisenberg model (Q). So 
we are left with the interaction between localized and itiner- 



the interaction term eventually reads, 



(13) 



Reintroduction of the lattice site dependence (Ri a ) we end up 
with the interaction term of the multiband-KLM: 



J 



'iam ^iam-a ^lama J 



(14) 
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FIG. 1: Supercell geometry for the EuO(lOO) film calculations, with 
the Europium atoms (large), the oxygen atoms (small) and the empty 
spheres (circles). The vertical lines indicate the surface of the EuO 
film. 



Here we have used the abbreviations 



Si x 4-7 7 ^ V 

'-'jam ~ b <Jj a 



(15) 



Compared to the conventional single-band KLM we have 
simply an additional band-summation in the exchange term. 
Note, however, the subbands are intercorrelated via the single- 
particle term Hd (§)■ 



B. Band structure calculations 



in the single-particle Hamilto- 



rnm 
ijafi 

has to contain the influences of all those interactions 



Thehopping integrals T f 
nian 

which are not directly covered by our model-Hamiltonian. For 
this reason we perform a band structure calculation accoijding 
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to the 'tight binding '-LMTO- AS A program of Anderson 
The resulting single-particle energies are taken as single- 
particle input for the partial Hamiltonian (|j). 

Test calculations for bulk-EuO have revealed the well 
known 'gap problem' which appears in connection with the 
strongly localized character of the 4f-'bands'. A 'normal' 
LDA calculation makes EuO a metal with the 4f-levels lying 
well within the conduction band. To overcome this situation 
we have investigated two-approaches for dealing with the lo- 
calized 4f moments (seeliil for the analogous case of Gd). In 
the first, the 4f moments are treated as localized core elec- 
trons, which naturally gives an insulator ground state, while 
the conduction band region consists predominaaJ-b of Eu-5d 
states. In the second, an LDA+U calculationESEl has been 
performed, which has the advantage compared to the former 
procedure, that the flat 4f bands explicitly appear right be- 
tween the 0-2p bands and the Eu-5d bands. In this approach, 
a meanfield approximation for the direct Coulomb U and the 
exchange interaction between the 4f electrons is added to the 
LDA-energy functional. The parameters U and J can be used 
to position the 4f levels correctly in the band gap between the 
0-2p and the Eu-5d levels. Besides the fact, that adjustable 
parameters U and J somehow corrupt the idea of a 'first prin- 
ciples' calculation, and since we consider the 4f-levels in our 
model study only as localized magnetic moments (Q), the 
much simpler LDA calculation has to be preferred in particu- 
lar when our study focuses on more complicated film geome- 
tries. For the band-structure calculations- for EuO films one 
has to employ a supercell geometry as depicted in Fig. [IJ As 
a result, the system has a super-layered structure with super- 
layers consisting of n consecutive EuO(lOO) layers followed 



by m layers of empty spheres. That means, we have a system 
of periodically stacked EuO n-layer films, isolated from each 
other by m layers of empty spheres. 

The number m of empty layers has to be chosen large 
enough to have truly isolated EuO films, but on the other hand 
also as small as possible to cut down on the numerical ef- 
fort. In Fig. ||, the 5d-partial density of states of EuO(lOO) 
monolayers (n=l) is displayed for different numbers m of 
spacer layers. The DOS of the EuO monolayer converges 
quite quickly as a function of m. For the following we have 
chosen m=5 which, according to Fig. ^, seems to be a reason- 
able approximation for the case of ideally isolated EuO films 
(rn — > oo). Fig.|3j shows the partial densities of states of the 
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FIG. 2: Partial densities of states for the Eu-5d bands of a EuO(lOO) 
monolayer calculated for different numbers m of spacer layers of 
empty spheres. 

5d bands of the center layer of films of various thicknesses and 
for bulk EuO. The respective curves converge as a function of 
film thickness. Comparing the partial Eu-5d DOS of the cen- 
ter layer (a = 10) of a 20-layer EuO(lOO) film with that of 
bulk EuO, we see that apart from an unimportant constant en- 
ergy shift the main features of the densities of states match. 
Consequently, the center layer at the 20-layer film can be re- 
garded as a bulk-like environment. The still existing slight dis- 
crepancy may be accounted for by the fact that the supercell 
geometry in the band structure calculations for the EuO films 
destroys the cubic symmetry present in bulk EuO. Therefore, 
the splitting of the d-states into t2 g and e g orbitals is not valid 
anymore. Furthermore, the center of a 20-layer EuO( 100) film 
of course can be only an approximation in bulk EuO. Thus we 
believe that the slight discrepancy between the partial Eu-5d 
densities of states of the center layer of the 20-layer film and 
of bulk EuO is acceptable. 

In Fig. ^ the partial Eu-5d densities of states of a 20-layer 
film are displayed for the two surface layers, a = 1,2, and 
for the two center layers, a — 9, 10. We see that the local 



FIG. 3: Partial densities of states for the Eu-5d bands of the cen- 
ter layers (a: layer index) for EuO films of various thicknesses n. 
Bottom: Bulk EuO 



density of states in the center of the film is constant, but dif- 
ferent from that in the vicinity of the surface of the film. The 
lower band edge of the Eu-5d bands of the surface layer com- 
pared to those of the center layer, which represent-a bulk-like 
situation, indicates the existence of surface statesO. This in- 
teresting point will be discussed a little bit later. When we take 
the results of the just-described band structure calculation as 
input for our finite temperature quasiparticle band structure 
calculation, then we are running into the well known 'double 
counting problem' of just the relevant interactions, once ex- 
plicitly in the model-Hamiltonian and besides that implicitly 
in the above band structure calculation. How we avoid this 
problem is explained in the next Section. 



HI. MANY-BODY EVALUATION 

After having fixed the model-Hamiltonian (Q) and its renor- 
malized single-particle energies in (|3|) via a full band structure 
calculation we have now to solve the many -body problem pro- 
voked by the Hamiltonian H. There are two partial problems 
to be solved, one concerning the magnetic properties of the 
localized 4f moments , the other dealing with the temperature- 
reaction of the conduction band states on the magnetic state of 
the moment system. Since we are interested in the description 
of the ferromagnetic semiconductor EuO, we can assume an 
empty conduction band (n=0), more strictly, a single electron 
in otherwise empty conduction bands. Therefore, the mag- 
netic order cannot be influenced by conduction electrons and 
the magnetic part can be solved separately. 
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FIG. 4: Partial densities of states for the Eu-5d bands for different 
layers a of a 20-layer EuO(lOO) film. 

A. Local moment system 

The system of localized /-moments is described by the ex- 
tended Heisenberg Hamiltonian (Q). Here we want to stress 
once more that the single-ion anisotropy constant Dq is small 
compared to the Heisenberg exchange interaction, Do <C 
J"f ■ By defining the magnon Green function 

P^(E) = ((st,Si,)) E , (16) 

we can calculate the /-spin correlation functions by evaluating 
the equation of motion 

EPf(E) = 2H 2 8 a0 {S* a ) 

+ {{[St a ,n f ]_;Sr^ E . (17) 

The evaluation of this equation of motion involves the de- 
coupling of the higher Green functions on its right-hand 
side, originating from the Heisenberg exchange term, and 
the anisotropy term using the Random Phasep&pproximation 
(RPA) and a decoupling proposed by LinesEJ, respectively. 
The details of the calculation can be found in a previous 
papeitj. For brevity we restrict ourselves to present here only 
the results. For theJayer-dependent magnetizations of the/- 
spin system we getO 

{S* a ) {l + <p a ? s+ \S-^ a ) + <pl s +\S + l + <p a ) 

h _ (1 + Va)2S+ l 

(18) 



6 



where 



N 



(k) 

e /3B 7 (k) _ 



(19) 



where, again, TV is the number of sites per layer and (3 = ir^f • 
The summation in Eq. (19) runs over the n poles E 7 (k) 



of the Green function P£ (E) and the Xaa-y(k) is the weight 
of the 7'th pole in the diagonal element of the Green function 
PQ a (E). The poles and the weights can be calculated from 
the RPA solution of Eq. (Jl7|): 



P^(E) = 2h 2 



with 



/ {Si) 







(S z n ) 



{El-A)-\ (20) 



^ = (D„$ Q + 2^ W{S*))5 afi - 2jf{S%). 

7 

(21) 

The $ Q come from the decoupling of the higher Green func- 
tion on the right-hand side of Eq. ( JlV[ ) and originates from ihe 
anisotropy term in the Hamiltonian H / according to LinesEl: 

2H 2 S(S + l)-3h(S z a )(l + 2<p a ) 



$ = 



(22) 



Having obtained the layer-dependent magnetizations (|l8|) 
and the ip a we can now express all the other /-spin correlation 
functions which we need for the following via the relations 

(S~S+) = 2h(S z a )^ a , (23a) 
{(S z a ) 2 ) = h 2 S(S+l)(S z a )(l + 2^ a ), (23b) 
((S z a f) = H 3 S(S+l) Va 

+h 2 (s z a )(s(s + i) + ^ a ) 

-h((S z a ) 2 )(l + 2<p a ), (23c) 
and the general spin-operator equality 



= k 2 s(s + 1) + z a hS z a (S z a ) 2 



(24) 



We shall see in the next Section that these/-spin correlation 
functions provide the whole temperature dependence of the 
electronic subsystem I 



B. Conduction bands 

Similarly as the magnetic part in the preceding Section the 
electronic part can be treated separately in the case of EuO. 
The reason is that magnon energies and exchange integrals, re- 
spectively, are smaller by some orders of magnitude than other 
energies in the system as the d-f exchange J or the conduction 
band width W. Accordingly, the respective terms (J? ■ , Dq) 
can be neglected when dealing with the electronic excitation 



spectra. That does not at all mean that there is no influence 
of the magnetic moments on the electronic quasiparticle spec- 
trum. That is done by the spin correlation functions of the last 
Section. 

Starting from the Hamiltonian of the electronic subsystem, 



n. — rid + ridf, 



(25) 



all physical relevant information of the system can be derived 
from the retarded single-electron Green function: 



Gram ( rp\ 



-tamer, ^jfim'a 



dte~t Et 



(0) 



■(26) 



Here and in what follows [.,.]+([.,.]_) is the anticommuta- 
tor (commutator). Conform to the two-dimensional transla- 
tional symmetry, we perform a Fourier transformation within 
the layers of the film, 

G^(E) = c ik(R *- R ^G™V(£), (27) 



where N is the number of sites per layer, k is an in-plane 
wavevector from the first 2D-Brillouin zone of the layers 
and Hi represents the in-plane part of the position vector, 
Ri a = Ri + r a . From Eq. (^7|) we get the spectral density 
by 



1 



S^JE) = - : ImG= a (E + i0+), 



(28) 



which is directly related to observable quantities within an- 
gle and spin resolved direct and inverse photoemission experi- 
ments. Finally, the wave-vector summation of S^J" (E) yields 
the layer-dependent (local) quasiparticle density of states: 

PZ(E) = -^J2 S ^(E). (29) 

k 

In the following discussion all results will be interpreted in 
terms of the spectral density ( ffi| ) and the local density of states 
©. 

For the solution of the many-body problem posed by Eq. 
(|25|) we write down the equation of motion of the single- 
electron Green function ( 26 ) 



Erpmva" frn 
ikc^y ^kj^ficr 



(30) 



The formal solution of Eq. ( j3C| ) can be found by introducing 
the self-energy M^ a {E), 



' C n 

' jpm'a 



5^ Mika\a(E)G™j7]3a(E)i 
k'ym" 



(31) 
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which contains all information about the correlations between 
the conduction band and localized moments. After combining 
Eqs. ( |3C| ) and ( |3~l| ) and performing a two-dimensional Fourier 
transform we see that the formal solution of Eq. ( ^o| ) is given 
by 



Gkcr(-E') = H (E I — Tk — M^jfi?)) -1 . 



(32) 



Here, I represents the (nM x nM) identity matrix, with 
M denoting the number of subbands of the conduction 
band system, and the matrices G\ (itT (E), Tk, and Mk„(E) 
have as elements the layer- and subband-dependent functions 
G^ a {E), Tg%', and M%$ a {E), respectively. 

To explicitly get the self-energy in Eq. ( plj ) we evaluate the 
Green function 



J 



' jpm'cr 



2 ^ cr iijotctpcr ' iijaotftcr J ' 

Here, the two higher Green functions, 



ikja-y(3cr \ 



T?m m { zp\ 
r ikja-yBcr\ nj ) 



^io> ^k-ym" —a i C 



f 



(33) 

(34) 
(35) 



originate form the two terms of the d-f Hamiltonian ( |l4| ) 
and will be referred to as the Ising and the Spin-flip 
function, respectively. Considering the equations of mo- 
tion for these two Green functions we encounter the two 



higher Green functions i[S? a Ck rfm ,"<T,'H<lt]-,Cjp m iJ)E and 

{([ S ia c k im ''-a,Hdf]-]c+ f3m , IJ ))E- Since we consider an 
empty conduction band the thermodynamic average in the 
Green functions has to be computed with the electron vacuum 
state | n = 0). From the definition of the t/-/Hamiltonian ( |l4] ) 
we then see that (n — 0| Hdf — and, accordingly, 



non-diagonal terms, (i, a, m) ^ (k, 7, m") and for the diag- 
onal terms, (i,a,m) — (k, 7,771"). In the first case we use 
a self-consistent so-called self-energy approach which results 
in a decoupling of the equations of motion. For the diagonal 
terms, (i,a) = (k, 7), this approach is replaced by a mo- 
ment technique which takes the local correlations better into 
account. 

The details of the approach have already been published 
in our previous paperll By introducing the multi-indices 

A, B, defined by (a, to) — ► A, (f3, m') — > B,..., 
the approximative solution forlhe electronic subsystem of the 
single-band KLM presented ir£l can be transferred to the case 
of the multi-band KLM, with the correspondence a = 4,j3 = 

B, — 

The reason for this straightforward transference is that the 
commutator relations which govern the solution inQ are invari- 
ant under the index transition^. g. 



[ciAa,H' t 



mbl 
df J- 



= ~ 2 ( Z <r S iA c iAo + <HA-a) (38) 



Thus there is no need to repeaLthe procedure once more. For 
details the reader is referred td/J 



IV. RESULTS FOR EUO-FILMS 

EuO crystallizes in the NaCl structure. Hence, the Eu 2+ 
ions are arranged on an fee lattice. For all the following 
film calculations, we assume the surface to be parallel to the 
fcc(100) crystal plane. 



[S ia , T~Cdf]_ c k-ym"o-'t c ~^f) m <a 



(( [ S iJi ftdf] _ Ck~fm"-<r]C 



+ 



0. 



Hence, for the equations of motion of the Ising and the Spin- 
flip function we get 



]T (E5 kl 6 



rpm m 
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(36) 
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(37) 



A. Magnetic properties 

The calculation of the magnetic properties of the EuO 
films has been done with the theory of Sect.III.A. We 
take the Heisenberg-exchange integrals from the experiment. 
Low temperature neutron scattering measuring the spin-wave 
dispersionE! -and being compared to the renormalized spin 
wave theorytl have revealed that nearest ( J\) and next-nearest 
neighbor ( J 2) exchange integrals have to be taken into ac- 
count: 



Ji_ 

k B 
•h 



0.625 K 
0.125 if 



(39) 



For an fee (100) film, we have therefore for the case of uni- 
form Ji and J2 within the film: 



On the right-hand side of these equations appear further higher 
Green functions which prevent a direct solution and require 
an approximative treatment. The treatment is different for the 
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FIG. 5: Dependence of the Curie temperature Tc on the single- 
ion anisotropy constant Do, calculated for bulk EuO. Dotted line 
represents the experimental value, Tc = 69, 33K. Inset: Mag- 
netization of bulk EuO as function of temperature, calculated for 
Do = 0.005/cb 



Here,Ai|i,Aij_ and A 2 ii denote,respectively, the positions of 
nearest neighbors within the same layer, within the adjacent 
layers and of next-nearest neighbors within the same layer one 
finds, 

2A M | = (1,1), (1,-1), (-1,1), (-1,-1) 
2Ai_l = (0,1), (0,-1), (1,0), (-1,0) 
A 2 |, = (0,1), (0,-1), (1,0), (-1,0) 
A 2± = (0,0) (41) 

It remains to fix in a reasonable manner the single-ion 
anisotropy constant Dq. For this purpose we have calculated 
the Curie temperature of bulk (!) EuO with J%, J 2 from ( fil| ) 
and variable Dq. The result is shown in Fig. || When p 2 - in- 
creases from to 0. IK Tc rises from about 66. 5K to 70K. Re- 
garding that J 1, J2 stem from a low-temperature fit, the agree- 
ment between the calculated Tc's and the experimental value 
of 69.33KU is remarkably good and that for all Do-values 
in the investigated region. We have chosen ^ = 0.05K 
for our calculations of film geometries. This is in reason- 
able agreement with the rather spare experimental values. 
According to ( p0[ ) and (|Tj) the model anisotropy energy is 
E a = —Do Q a with the temperature-dependent anisotropy 
coefficient <£> Q from Eq. (|22|). It can be seen that E a mono- 
tonically decreases as function of (S*) to E a = — 6D0 at 
{S z a ) = S. That means at T = we have w -0.3 which 
is the same order of magnitude as the experimental value in 
refH 

Fig. H shows the temperature and layer-dependent magne- 
tizations of EuO(100) films for thicknesses from n = 1 to 
n = 20. For all temperatures and for all film thicknesses the 
(Sa) increase from the lowest magnetization of the surfaces 
monotonously towards the highest magnetizations in the film 
centers. That can qualitatively be explained by the lower coor- 
dination number of the surface atoms, Zsj cc (ioo) = 4, com- 
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T/K 
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FIG. 6: Layer-dependent magnetizations, (S„), of EuO(100) films 
as functions of temperature for various thicknesses n. (S Q ) always 
increases monotonously from the surface layer towards the film cen- 
ter. Inset: Curie temperature as a function of film thickness. 



pared to the bulk coordination number, Zh./ cc = 12. The inset 
of Fig. H shows the dependence of the Curie temperature on 
the film thickness n compared to the bulk-Tc calculated with 
the same parameters J\, J2, and Do (see inset of Fig. |J). Start- 
ing from a rather low value for the monolayer Tq steadily in- 
creases with n reaching the bulk-value; for n > 20. The curve 
strongly resembles that of Gd fiimstij. Corresponding mea- 
surements for EuO are unknown to us. 

All the higher spin correlation functions (S^S^), ((S^) 2 }, 
((S^) 3 ), . . . follow from yia (|l|), @, (jSj) . Exam- 
ples are plotted in Fig. || of ref.Q. Together with (S^) they 
provide the necessary temperature information needed to cal- 
culate the temperature-dependent electronic structures of EuO 
films, which are presented in the next Section. 



B. Temperature-dependent electronic structure 

In order to get the temperature-dependent band structures 
of the EuO(100) films we have to combine the LSDA calcu- 
lations, described in Sect.II.B, with the many-body evaluation 
of the multiband-Kondo hffXce model (Sect.III.B). It has been 
shown in previous worksatHJ that the d-f exchange interaction 
of the KLM can be treated exactly for T=0 in the case of empty 
conduction bands. While the |-spectrum turns out to be rather 
complicated exhibiting interesting correlation effects, the f- 
spectrum is simple being only slightly shifted towards lower 
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energies by a constant energy amount of i J S. So we can di- 
rectly use the TB-LMTO obtained k-dependent hopping ma- 
trices for the j-electron as input for the single-particle Hamil- 
tonian H c i (Q). The problem of double-counting of relevant 
interactions, here the d-f-exchange interaction, which usually 
occurs when combining first-principles and model calcula- 
tions, is then elegantly avoided. The temperature comes into 
play via the layer- and temperature-dependent f-spin correla- 
tion functions, discussed in the last section, which appear in 
the many -body evaluation of the KLM (Sect. )TIB ). To avoid 

ambiguities, we use in the following as temperature-parameter 

(s z ) 

the reduced magnetization of the center layer of the re- 
spective EuO (100) film. The same temperature may produce 
in the other layers different magnetizations. 

Table | gives some examples how the different values of the 
magnetization are related to the respective temperatures 
for EuO (100) films. The only parameter of our theory, the. 
exchange coupling J, is chosen from bulk-EuO calculationsEl 
to yield the experimentally observed 'red shift' of the lower 
conduction band edge upon cooling down from T — Tq to 
T = 0: 



J = 0.25 eV 



(42) 




(S,)/S 


bulk 
EuO 


EuO(100) films (number of layers) 


1 
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10 


20 
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. 75 


45.48 


10 . 52 


21.56 


37 . 41 


43 .23 


45.14 


. 5 


59.76 


13 . 95 


28.51 


47 .51 


55 .59 


58 . 92 


0.25 


66.71 


15 .56 


31.82 


52 .02 


61 . 02 


65.04 


(To) 


68 .84 


16.06 


32.83 


53.36 


62.70 


66.73 



TABLE I: Temperatures at which the 4f magnetization of bulk EuO 
and of the center layers of EuO (100) films has assumed the values 
in the left column (see Figs, pi 




FIG. 7: Spin-dependent spectral densities of the Eu-5d bands of a 
EuO(100)monolayer for different magnetizations (%-} (see Table 



Energy (eV) 

FIG. 8: Temperature-dependent densities of states of the Eu-5d bands 
of a EuO(100) monolayer for different values of the 4f magnetization 
(^) (see Table!). 



Let us start the discussion with the spectral densities and the 
quasiparticle densities of states (Q-DOS) of the Eu-5d bands 
for a EuO(100) monolayer (n=l) at varying magnetizations 
of the 4f moment system. Fig. [7] shows the quasiparticle 
band structure as a density plot. The degree of blackening is a 
measure of the peak-height of the respective spectral density. 

Broad, washed out lines mean quasiparticles with relatively 

( s z ) 

low life-times. In ferromagnetic saturation ( 1 g- i = 1) the j- 
spectrum consists of stable quasiparticles. As explained above 
this is the exactly solvable limiting case. The spectrum is iden- 
tical to the LSDA input only slightly shifted by — |j 5". The 
|-spectrum is also exactly solvable exhibiting, however, clear 
correlation effects as, e.g., a splitting of the dispersion around 
4eV near the -M-points. As discussed in detail for a model 
system in refBLl the two excitations belong to two different 
microscopic interaction processes. The high-energy part is 
due to the formation of a magnetic polaron (propagating elec- 
tron 'dressed' by a virtual cloud of magnons), the low-energy 
branch is due to scattering processes (magnon-emission by the 
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V 



finite temperatures (demagnetizations) magnons are available, 
absorption processes lead to damping of the quasiparticle dis- 
persion in the |-spectrum, too. The overall exchange splitting 
reduces with increasing temperatures and the two spin spec- 
tra are becoming more and more similar. Finally, in the limit 
of {S z ) — ► 0, T — ► T c the lack of any 4f-magnetization 
removes the induced spin-asymmetry in the 5d bands (Fig. f7j). 

The temperature-dependent Q-DOS of the EuO-5d bands 
in a monolayer exhibits a strong induced exchange splitting 
at T=0 (Fig. ^|), where the increase of temperature from T=0 
to T = Tq leads to a convergence of the two spin parts. It 
can clearly be seen that the temperature dependence of the 
exchange-split spectra goes beyond the Stoner picture, where 
the bands move rigidly towards each other when the temper- 
ature is increased. The reason lies predominantly in the spin- 
flip processes between 5d- and 4f- system. Such correlation 
effects take care for the non-rigid shift. By comparison with 
the f-Q-DOS for ferromagnetic saturation ((S z ) — S), the 
shape of which is identical to the LDA-renormalized single- 
particle input, one recognizes that in the paramagnetic phase, 
too, correlation effects remarkably alter the spectrum. 

Fig. U displays the spectral density of a EuO double layer 
for three different magnetizations of the 4f-moment system. 
The same which has been said above about the temperature- 
dependent tendencies of the monolayer EuO(lOO) spectra ap- 
plies to the spectra of the double-layer in Fig. g. However, 
due to the higher number of bands which originates from the 
hybridization of the two EuO(lOO) layers the temperature- 
dependent evolution of the bands becomes less clear for the 



FIG. 9: Spin-dependent spectral densities of the Eu-5d bands of an 
EuO(lOO) double layer (n=0) for different values of the 4f magneti- 
zation (%-} (see Table |). 



!-electron). The magnetic polaron can be considered as due 
to a polarization of the neighboring 4f spins by the conduction 
electron via d-f exchange. Under certain parameter constella- 
tions this may even leadJg a bound state,e.g. to a quasiparti- 
cle with infinite lifetimeBfl. This exchange-induced splitting, 
which can also be observed in other parts of the T=0-|-band 
structure, is a typical correlation effect, which by no means 
can be accounted for in an LSDA calculation. 

An exchange splitting of the full |- and j spectra against 
one another is clearly visible. Note that these effects appear 
although the bands are empty, more strictly there is one elec- 
tron ('test electron') in the otherwise empty 5d-bands. They 
are therefore certainly induced by the 4f system and mediated 
by the d-f exchange interaction. 

For finite temperature, intermediate magnetizations, -^U- = 
0.5, the f-dispersion, too, starts to wash out indicating spin- 
flip processes. Magnon emission by a ^-electron is completely 
equivalent to magnon absorption by a ^-electron, with one ex- 
ception: magnon absorption is possible only if there are any 
in the spin system. This is not the case in ferromagnetic sat- 
uration {(S z ) = S). That is the reason why at T=0K the f- 
spectrum appears so much simpler than the j-spectrum. For 



<S Z >/S = 1, 0.75, 0.5, 0.25, 




-2 2 4 6 8 

Energy (eV) 



FIG. 10: Local densities of states of the Eu-5d bands of the first, sec- 
ond and center layer of a 5-layer EuO(100) film for different values 
of the center layer 4f magnetization (see Table |). 



11 



<S >/S = 1, 0.75, 0.5, 0.25, 




Energy (eV) 



FIG. 1 1 : The same as in Fig. [u] but for a 10-layer film 



EuO(lOO) double layer. 

Figs. [l(] and |ll| show the temperature-dependent Q-DOS 
of, respectively, a 5-layer and a 10-layer EuO(lOO) film. Dif- 
ferent from the case of the monolayer and the double layer, for 
these films the densities of states depend on the layer index a. 
Furthermore, the magnetization of the 4f-moment system is a 

layer-dependent entity (see Fig. n). Here, the temperature pa- 

(s z ) 

rameter ■^- L refers to the magnetization of the center layer of 
the film. The respective temperatures are given in Table [| 

With increasing film thickness, the Q-DOS of the center 
layer becomes more and more similar to the Q-DOS of bulk 
EuO. As a result, the centers of gravity of the surface layer 
bands, on the one hand, and of the center layers, on the other 
hand, move away from each other, with the latter positioning 
at higher energies. This tendency leads to the existence of 
EuO(lOO) surface, states, which are discussed in detail in an 
additional papertl 

There it is shown that for T = Tq the surface at the F- 
point lies some 0.8 eV below the bulk conduction band edge. 
This distance does not change with temperature, t.e. the f part 
of the surface state exhibits the same red shift of about 0.3 
eV as the 5d conduction band edge when cooling down to 
Tcrr= OK. Since the paramagnetic 4f-5d gap amounts to 1.12 
eVtl it can be speculatedo that the f-surface state closes the 
gap at very low temperatures. Giving rise to a temperature- 
induced surface insulator-halfmetal transition. The Eu-ion 
then remains in an intermediate valence fluctuating between 
the magnetic 4f 7 -state (J = S = |) and the non-magnetic 
4f 6 -state (J = 0). Interesting dynamic effects are to be ex- 
pected since the non-magnetic 4f 6 -state automatically cancels 
the red shift pushing the |-surface states by 0.3 eV to higher 



energies. Electrons which will reenter the 4f level creating 
again the magnetic 4f 7 configuration, and so on. 



V. CONCLUSIONS 

By combination of a 'first-principles' band structure cal- 
culation within the frame of density functional theory with a 
many-body evaluation of a proper theoretical model we have 
derived the temperature-dependent electronic structure and 
the magnetic properties of a real 4f system (EuO) with film 
geometry. The theory provides the properties of films of vary- 
ing thicknesses. The film geometry provokes a superposition 
of electron correlation effects with the sheer geometrical de- 
pendence of the spectra on film thickness and layer index. 

'Local-moment' systems like EuO get their striking tem- 
perature dependencies and correlation effects by an exchange 
interaction between localized magnetic moments, due to 4f 
electrons of the Rare Earth ion, and quasi-free (5d,6s) conduc- 
tion electrons. The situation is rather convincingly modeled 
by the Kondo-lattice (s-f, s-d) model (KLM). We proposed a 
way how to generalize the original single-band KLM to a re- 
alistic multi-band model. The respective many-body problem 
is approximately solved for the electronic subsystem mak- 
ing use of a moment-conserving decoupling approximation 
(MCDA) for suitably defined Green functions. The method 
gains substance by the fact that the finite-temperature theory 
evolves continuously from the exactly solvable and nontrivial 
T=0 case of ferromagnetic saturation. 

The mentioned T=0-limiting case is also used to connect 
the multiband-model treatment with a TB-LMTO-ASA band 
structure calculation for thin EuO(100) films. Circumvent- 
ing the much discussed double-counting problem the spin-up 
Eu-5d bands have been taken as input for the hopping matrix 
of the multi-band-KLM Hamiltonian. The reason is that the 
exact T=0-KLM solution shows for ferromagnetic saturation 
that the f -quasiparticle spectrum is only rigidly shifted com- 
pared to the 'free' dispersions. A double counting of the d-f 
exchange interaction is therefore excluded. 

The film geometries are accounted for by a supercell con- 
struction, where consecutive EuO(100) films are isolated from 
each other by a stack of spacer layers. The connection to the 
model calculation employs a decomposition into subbands, 
that respects the symmetry of the Eu-5d orbitals. 

Already for the case of ferromagnetic saturation of the 
4f-moment system, the spin-down spectra of the EuO films 
exhibit strong deviations from the single-electron LSDA re- 
sults demonstrating strong correlation effects. These corre- 
lation effects in the j-spectra provoke a strong temperature- 
dependence. The |-spectra, which at T=0 correspond to the 
t-LSDA band structure results, develop similar correlation ef- 
fects for finite temperatures. The loss of 4f magnetization for 
T — ► Tq leads to a convergence of the spectra for both spin 
directions. 

The semiconductor EuO fits the physics of the KLM for 
the special case of an empty conduction band. The extension 
to non-zero band occupation allows to study the prototypi- 
cal 'local moment' metal Gd. Concerning the aspect of re- 
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